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A PROGRAM FOR CALCULATING EXPANSION -TUBE 
FLOW QUANTITIES FOR REAL -GAS MIXTURES AND COMPARISON 
WITH EXPERIMENTAL RESULTS 

By Charles G. Miller HI 
Langley Research Center 

SUMMARY 

A computer program written in FORTRAN IV language is presented which deter- 
mines expansion -tube flow quantities for real test gases CO 2 , N 2 , O 2 , Ar, He, and H 2 , or 
mixtures of these gases, in thermochemical equilibrium. The effects of dissociation and 
first and second ionization are included. Flow quantities behind the incident shock into 
the quiescent test gas are determined from the pressure and temperature of the quiescent 
test gas in conjunction with (1) incident -shock velocity, (2) static pressure immediately 
behind the incident shock, or (3) pressure and temperature of the driver gas (imperfect 
hydrogen or helium). The effect of the possible existence of a shock reflection at the 
secondary diaphragm of the expansion tube is included. Expansion -tube test-section flow 
conditions are obtained by performing an isentropic unsteady expansion from the condi- 
tions behind the incident shock or reflected shock to either the test -region velocity or the 
static pressure. Both a thermochemical -equilibrium expansion and a frozen expansion 
are included. Flow conditions immediately behind the bow shock of a model positioned at 
the test section are also determined. Results from the program are compared with pre- 
liminary experimental data obtained in the Langley 6 -inch expansion tube. 

A listing of the program is presented along with a description of required inputs, a 
flow chart, and sample data printouts. Some typical solutions are presented for CO 2 , Ar, 
and air and for Mars, Venus, and Jovian atmospheric models. 

INTRODUCTION 

Several experimental studies to be performed in the Langley 6 -inch expansion tube 
will use various gas mixtures as test media at hypersonic conditions. Before such studies 
are performed, it is essential to ascertain the theoretical performance of the expansion 
tube for the gas mixture being tested. The wide range of flow conditions and the very 
short test times of the expansion tube impose rather stringent requirements on facility 
instrumentation. Thus, in preparing the facility for a test, it is necessary that the magni- 
tude of the physical quantities to be measured be known to within reasonable limits. 



Following a test, it is desirable to have convenient means for determining expansion -tube 
flow quantities from the measured flow quantities. 

Relatively little effort has been directed towards the calculation of expansion -tube 
flow quantities for arbitrary gas mixtures. One such study, reported in reference 1, is 
basically a shock-tube program which can generate expansion -tube flow conditions. How- 
ever, several effects which may markedly alter calculated flow quantities are not included 
in the program of reference 1. Also, the program of reference 1 requires inputs other 
than the actual flow quantities measured during a test. 

A primary purpose of the present study is to fulfill the need for a convenient, versa- 
tile program for determining expansion-tube flow quantities in real gases CC>2, N£, C>2, 

Ar, He, and H2 or mixtures of these gases from several combinations of routinely mea- 
sured flow quantities. Included in this program are the effect (ref. 2) of an imperfect 
(intermolecular force) driver gas (helium or hydrogen), the possible existence of a shock 
reflection at the secondary diaphragm (ref. 3), and the calculation of both a thermochemical - 
equilibrium expansion and a frozen expansion. Flow quantities immediately behind the 
bow shock of a test model, including stagnation -point heat -transfer rate, are calculated. 

SYMBOLS 

The International System of Units (SI) is used for all physical quantities in the pres- 
ent study. Conversion factors relating SI Units to U.S. Customary Units are given in 
reference 4. 

a speed of sound, m/sec 

B second virial coefficient, m^/kg 

2 

C third virial coefficient, (m^/kg) 

CpWo/R nondimensional specific heat at constant pressure 

c v W 0 /R nondimensional specific heat at constant volume 

h specific enthalpy, m^/sec^ (j/kg) 

h Q standard heat of formation at T = 0 K, j/kmol 

M Mach number, U/a 
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M s 

%e 

P 

q 

R 


u 

u< 

w 

w, 

Xi 

Yi 


o 


E 


incident -shock Mach number, Ug^/a 

Reynolds number per meter, pU/p 

pressure, N/m^ 

heat -transfer rate, W/m^ 

universal gas constant, 8.31434 kJ/kmol-K 


r nose radius, m 

sW 0 /r nondimensional specific entropy 

T temperature, K 


velocity, m/sec 

velocity of incident shock, m/sec 
molecular weight, kg/kmol 

molecular weight of undissociated gas or gas mixture, kg/kmol 

mole fraction, kmol of species i per kmol of mixture 

number of kmoles of species i per mass of mixture, kmol of species i 
per kg of mixture 

compressibility factor, pW 0 /pRT 

number of kmoles of dissociated gas mixture per number of kmoles of undis- 
sociated gas mixture, W Q/ /W 

ratio of specific heats 


Vi 


isentropic exponent, fjLi2I_R] 

\a log p/sw o /R 

parameter defined in equation (25), kg/^i^/^-sec^N/m^)'*^ 
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(i viscosity, N-sec/m^ 

p density, kg/m^ 


Subscripts: 


A 

act 

f 

1 

mix 

perf 

ref 

w 

1 

2 

2, s 

3 

4 

5 

5,s 

5jt 


denotes region (IT) for no standing shock at secondary diaphragm or 
region (2^s) for a standing shock 

active or available 

frozen 

individual species 
mixture 

perfect -gas behavior 
reference state 
model -wall conditions 

state of quiescent test gas in front of incident normal shock 

state of test gas behind incident normal shock in intermediate section 

state of test gas behind shock at secondary diaphragm 

state of expanded driver gas 

driver -gas conditions at time of primary -diaphragm rupture 

state of test gas flow in acceleration section (free -stream conditions) 

static conditions behind bow shock of model positioned at expansion -tube 
test section 

stagnation conditions behind bow shock of model positioned at expansion -tube 
test section 
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10 


state of quiescent acceleration gas in front of incident normal shock 


20 state of acceleration gas behind incident normal shock in acceleration section 

Approximate value is denoted by superscript 

ANALYSIS AND PROCEDURE 

Before the procedures for determining expansion -tube flow quantities are discussed, 
a brief description of the Langley 6 -inch expansion tube and its operating sequence is 
given. This facility is basically a cylindrical tube, having a 15.24-cm inside diameter, 
divided by two diaphragms into three sections (fig. 1(a)). The most -upstream section is 
the driver or high-pressure section. This section is pressurized at ambient temperature 
with a gas having a high speed of sound, such as hydrogen or helium. (Greater operation 
efficiency is realized with gases having high speed of sound.) The pressure and speed of 
sound of the driver gas may be increased further by heating the gas with a 3-MW resis- 
tance heater or by utilizing an arc discharge into the gas from a 10-MJ capacitor bank. 

The intermediate section is usually referred to as the driven section. This section is 
evacuated and then filled with the test gas at ambient temperature. The most -downstream 
section is denoted as the acceleration or expansion section. This section is also evacuated 
and is generally filled with helium at a low pressure and at ambient temperature. For 
unheated or resistance heating of the driver gas, the driver and driven sections are sepa- 
rated by a double -diaphragm apparatus capable of withstanding a maximum pressure dif- 
ferential of 68.95 MN/m^. (By controlling the pressure level in the small chamber 
between these diaphragms, the time of diaphragm rupture can be controlled.) For arc 
heating, a single diaphragm is used between the driver and driven sections. A weak, low- 
pressure diaphragm (secondary diaphragm) separates the driven and acceleration sections. 
The test section and the model are located at the downstream exit of the acceleration 
section. 

The operating sequence, which is shown schematically in figure 1(a), begins with the 
rupture of the primary or high-pressure diaphragm. A primary shock wave propagates 
into the static test gas and an expansion wave propagates into the driver gas. The shock 
wave then encounters and ruptures the secondary or low-pressure diaphragm. The flow 
energy lost in rupturing this diaphragm results in an upstream -facing shock wave reflected 
from the diaphragm (ref. 3). This is shown schematically in figure 1(b). A secondary 
shock wave propagates into the low-pressure acceleration gas while an upstream - 
expansion wave moves into the test gas. In passing through this upstream -expansion wave, 
which is being washed downstream since the shock-heated test gas is supersonic, the test 
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gas undergoes an isentropic unsteady expansion resulting in an increase in the flow veloc- 
ity (ref. 5 ). Testing takes place in the flow that has passed through the expansion. This 
test region is denoted as region (If) in figure 1. 

Determining Flow Quantities in Region (If) 

The first step in determining expansion -tube flow conditions is to calculate the flow 
quantities in region (If) (fig. 1 ). In the Langley 6 -inch expansion tube, the quiescent -test - 
gas pressure p^ and temperature Tj (ambient) are known, as is the test -gas composi- 
tion Xj. The incident -shock velocity U s i is routinely inferred from microwave mea- 
surements and from the response of instrumentation spaced at known intervals along the 
driven section, and the pressure behind the incident normal shock in the intermediate 
section P2 is measured directly.. The flow quantities in region (ff) may be determined 
by using either Ugj or p2 as an input, in conjunction with p^, Tj, and Xj. Also 
considered is the case where p^ and T4 are known for the hydrogen or helium driver 
gas and are used in conjunction with pj, Tj, and Xj. 

The basic equations required in determining flow quantities in region (T) are the 
conservation relations for mass, momentum, and energy for a normal shock moving 
through region (T), which are 


Pl U s,l - P 2 ( U s,l " U 2) 

(1) 

P 1 +P 1 U S,1 2 " p 2 +p 2( u s,l - u 2) 2 

(2) 

hj + | u s,l 2 = h 2 + i(U s>1 - U 2 ) 2 

(3) 


and the equation of state (that is, source of thermodynamic properties for real-gas mix- 
tures). Thermodynamic properties for CO2, N2, O2, Ar, He, and H2 gases or mixtures 
of these gases in thermochemical equilibrium are obtained by the procedure of refer- 
ences 6 and 7 . The procedure of references 6 and 7 is included in the present study as a 
subroutine and is referred to herein as ROGO. An iterative -interpolation scheme which 
enables determination of thermodynamic properties from ROGO, for input combinations 
of h, p, sW 0 /R, and p, is also a subroutine and is referred to herein as FIND (I). (See 
appendix A for a detailed discussion of the source of thermodynamic and transport prop- 
erties for the arbitrary gas mixtures used in the present study.) 

Since p-^ and Tj are known, the corresponding thermodynamic quantities in 
region (T) (that is, p^ and hj) appearing on the left sides of equations ( 1 ) to ( 3 ) are 
obtained from the perfect -gas relations 
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( 4 ) 


_PlW 0 

P 1 “ RTi 

and 


V\ R VXj ih 0 

1 = — — r — — T-. + ) 1 > 1 

1 - 1 Wo 1 4 w 


The corresponding is obtained from the relation 



where 



(5) 


( 6 ) 


(7) 


For the Langley 6-inch expansion tube, Tj is the ambient temperature and pj is gen- 
erally less than 1.0 MN/m 2 . For these conditions, imperfect -gas effects are negligible. 
The iterative procedures for determining the quantities in region (2) from equations (1) 
to (3) for inputs U g j, p 2 , and p^ and T4 are now discussed individually. For all 
three cases, the flow in region (IT) is assumed to be in thermochemical equilibrium. 

Case where U g ± is known.- When U g ^ is known, the terms appearing on the 

left sides of equations (1) to (3) are known. An initial estimate of p 2 = lOp^ is made 
and the corresponding values of U2, p 2 , and h 2 are obtained. The p 2 and h£ are 
used as inputs to FIND (2) to obtain a value of p 2 . This p 2 from FIND (2) is compared 
with the initial guess of p 2 and, if it is not within 0.1 percent, the p 2 from FIND (2) is 
used in equations (1) to (3) to obtain upgraded values of U2, P 2 , and h£. This iterative 
procedure, commonly referred to as the method of successive approximations, is con- 
tinued until successive values of p 2 are within 0.1 percent. 

Case where p 2 is known.- When p 2 is known, an initial estimate of p 2 is made. 
The corresponding U s 1 is obtained from equations (1) and (2) in the form 


U 


s,l~ 


P 2 - Pi 


-ll/2 



( 8 ) 


7 



U2 is found from equation (1), and h2 is found from equation (3). The P2 and h2 
are used in FIND (2) to obtain a corresponding value of pg. The method of successive 
approximations is employed on pg until successive values are within 0.1 percent. 

Case where P4 and T4 are known. - The thermodynamic properties in region (?) 


are determined from the imperfect -gas relations of appendix B. The unsteady expansion 
which occurs from region (?) to region (3) (fig. 1(a)) upon rupture of the primary dia- 
phragm is assumed to be isentropic. An array of thermodynamic properties p, h, Z, 
a, and p is generated in region (3) by varying T3 from a maximum value (T4) to a 
minimum value and by knowing S3W3/R (S3W3/R = S4W4/R). The U3 corresponding 
to each value of T3 is obtained by using the differential equation for a one -dimensional, 
unsteady expansion. Across an upstream -facing unsteady -expansion wave, the velocity 
increment is related to the thermodynamic properties by the differential equation (ref. 5) 



Integrating equation (9) between regions (?) and (?) gives 



( 10 ) 


where U4 is zero. The U3 is found by performing the integration of equation (10) 
numerically. Thus, a unique correspondence between U3 and p^ at constant s^W 
is obtained. 


4 


/ R 


By varying U s> i over a range, an array of U2 and P2 may also be obtained. A 
solution is found by requiring that Pg equal Pg and U3 equal U2, as is illustrated 
in the sketch 


U 2 


and 



The solution is the intersection of the U2,P 2 curve and the Ug,pg curve. The other 
thermodynamic properties in region (?) are obtained from ROGO, where T 2 is obtained 
by interpolation and U Sj i is determined from equation (1). 
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Determining Flow Quantities in Region (£^s) 

As discussed in references 3 and 8, the flow energy lost in the rupture of the sec- 
ondary diaphragm must result in an up stream -facing shock wave reflected from the 
diaphragm. In reference 9, this shock-wave reflection is observed to be, for a short 
period of time, approximately a standing wave, as illustrated in figure 1(b). Therefore, 
the possible existence of a standing normal shock at the secondary diaphragm, due to 
shock-wave reflection from this diaphragm, is considered. The conservation relations 
for this standing shock are 

( 11 ) 


and 


p 2 u 2 - P2 >s U 2,s 

P2 + P 2 U 2 2 = p 2,s + p 2,s U 2,s 2 


, In 2 , 1 1T 2. 

h 2 + 2 U 2 _h 2,s + 2 U 2 ,s 


( 12 ) 


(13) 


Since the quantities in region (IT) are known (that is, have been calculated previously), 
equations (11) to (13) are solved by the method of successive approximations to yield the 
quantities behind the standing normal shock (that is, flow quantities in region (2^s)^. 

As in region (^), the flow in region (2^s) is assumed to be in thermochemical equilibrium. 


Determining Flow Quantities in Region (IT) 

Upon rupture of the secondary diaphragm, an upstream -expansion wave moves into 
the test gas. In passing through this upstream -expansion wave, the test gas undergoes an 
isentropic, unsteady expansion to region (IT). Integrating equation (9) between regions (A) 
and (jf) (where region (A) denotes region (If) for no standing shock at the secondary 

diaphragm or region (2^s) for a standing shock and where the conditions are as yet 

gives 


unknown in region f 5 

au = u 5 - u A = -j 


*5 /dll' 


h A 


V a /s a W 0 /R 


(14) 


Since p^ and U5 are measured in the Langley 6-inch expansion tube, both quantities 
are considered, individually, as inputs necessary for the solution of equation (14). As is 
typical of high -enthalpy facilities, the assumption of thermochemical -equilibrium flow is 
subject to question. Hence, both equilibrium -flow and frozen-flow cases are considered 
in the present program. 


Thermochemical-equilibrium unsteady expansion . - If the case where U5 is an 
input is considered first, then AU of equation (14) is known. The expansion in velocity 
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from region (A) to region (?) corresponds to a decrease in pressure p between these 
regions. Thus, with s^W 0/ /R and a range of p from an upper limit of p^ to a lower 
limit of Pg j (calculation of Pg ^ is discussed subsequently), FIND (3) is used to gen- 
erate a plot of 1/a against h for the real-gas mixture. The temperature limits, as 
required by FIND (I), used in generating this plot are T^ as an upper limit and Tg^ 
as a lower limit. Equation (14) is integrated numerically (Simpson’s rule) between the 
known limit h^ and the unknown limit hg. The hg is the value of h which equates 
the integral of equation (14) to AU. The corresponding Pg and Tg are obtained by 
interpolation and the other thermodynamic quantities in region (?) are obtained from 
ROGO. 


For the case of an equilibrium expansion where Pg is an input, s^W 0 y/R and Pg 
are used as inputs in FIND (3) to obtain the corresponding thermodynamic properties in 
region (?). With the limits of integration known, the integral of equation (14) is evaluated 
numerically. Once a value for AU is obtained and with known, U 5 may be found. 

Additional free-stream( region (?)) quantities that are calculated are Mg = Ug/ag and 

N Re,5 = p 5 U 5/^5* 

Frozen unsteady expansion .- Frozen flow is defined herein as that in which the 
vibrational energy and chemistry remain unchanged during the expansion process. In the 
present program, this freezing of the vibrational energy and chemistry is assumed to 
occur in region (A) . Hence the energy in region (A) may be viewed as consisting of an 
active or available part, which provides the energy for flow expansion, and a frozen or 
nonavailable part. In the frozen expansion, y^ remains constant and the test-gas mix- 
ture behaves as a perfect gas. By using the perfect -gas relation 


dh 



a da 


equation (14) may be evaluated in closed form to yield 


(15) 


U 5,f - U A - 5^T(*A,f - *5,f) < 16 > 

By freezing (that is, maintaining constant) the vibrational energy, CpWq/k becomes 7/2 
for a diatomic molecule and 9/2 for a triatomic molecule. The CpW^R for an atom 
is 5/2. By using the perfect-gas relation 


c v W o _ c p W o 


R 


R 


( 17 ) 


the specific-heat ratio y^ is obtained from equation (7) as 


© 
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The speed of sound in region (A) is determined from the relation 

/ R 7 * T \ 1//2 
a A,f " ( y f W 0 Z A T Aj 

and the active enthalpy is determined from the relation 


( 18 ) 


. a A,f , i 


I X i,A h o,i 


A, act ' - 1 


W r 


(19) 


Ionized species in region (A) are not included in the calculation of h A act from equa- 
tion (19). The values of h 0 for species N, N 2 , O, 0 2 , NO, C, CO, C0 2 , and Ar are 
given in reference 6 and those for He, H, and H 2 are obtained from reference 10. The 
enthalpy is given by the relation 


h A,f = h A ' h A,act ^ 

For the case where U5 (that is, U 5 ^ is known, a & f may be found from equation (16) 
and the corresponding frozen-flow thermodynamic quantities in region (J) are deter- 
mined from the isentropic perfect -gas relations of reference 11. For the case where p 5 
(that is, p 5 A is known, a 5 f is determined from the isentropic perfect -gas relation 


l 5,f 


= a 



l/2y f 


(21) 


Corresponding frozen-flow quantities in region ( If) are determined similarly, and with 
a^ | known, Ug f may be obtained from equation (16). 


Determining Flow Quantities Behind Bow Shock of Test Model 

In most tests in the Langley 6 -inch expansion tube, a test model (for example, a 
pitot probe) is positioned at the test section (tube exit). Hence, it is desirable to deter- 
mine the flow quantities behind the bow shock in front of a blunt test model. In the pres- 
ent study, flow quantities are calculated for the static conditions immediately behind the 
portion of the bow shock that is essentially a standing normal shock (region (5^s)j and 
for the stagnation conditions (point (5^t)j. These conditions are illustrated in the follow- 
ing sketch: 
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The conservation relations for a standing normal shock about a blunt body in the 
expansion -tube test section are given in equations (11) to (13) where the subscripts 2 
and 2,s are now replaced by 5 and 5,s. In the present study three combinations 
concerning the flow chemistry in the free-stream region (region (lT)j and the post -normal- 

shock region (region (5^s) and point (sjt)) are considered. These combinations are 
designated cases (1) to (3) as follows: 


Case 

Region (b) 

Region (5^s) and point 

© 

n 

Equilibrium 

Equilibrium 


■ 

Frozen 

Equilibrium 


mm 

Frozen 

Frozen 



For case (1) or (2), where post-normal-shock flow conditions are assumed to be in ther- 
mochemical equilibrium and free-stream conditions are either in thermochemical equilib- 
rium or frozen, the conservation relations are solved by iteration for the static conditions 
immediately behind the shock. In this solution, the flow quantities in region (IT) are con- 
sidered known and the initial estimate of p^ is taken to be 10 p^. For the post -normal- 
shock stagnation conditions, solutions are found by assuming that the flow region from 
immediately behind the normal shock to the stagnation point is isentropic (that is, 


s 5,t W ' 


o/R ■ s 5 ,s W o/ R > 


Now, the variations in conditions from immediately downstream 
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and 


of the shock to the stagnation point are relatively small (that is, 5 s ~ 5 t 

z *5 s * z *5 t)- Thus > the g as in this region may be considered to behave as a perfect 
gas. The p 5 t can be accurately determined from the perfect -gas isentropic relation 
(ref. 11) 


_ p + r E.5,s 

>,t p 5 f sy 2 


-1 


M 


y E,5,s 

i y B,5,s- 1 


5,s 


Now, for case (1), hg t is determined from the energy relation (eq. (13)) to be 


h 5,t = h 5 + ! U 5 2 


whereas, for case (2), the energy equation becomes 

h 5,t = h 5,s + | u 5,s 2 a h 5,f + h A,f + | u 5,f 2 


( 22 ) 


(23) 


(24) 


that is, for case (2), the nonavailable energy in the expansion process or free stream is 
made available to the post-normal-shock flow. The p 5 t from equation (22) and the 
h 5 t from equation (23) or (24) are used as inputs to FIND (2) to obtain the corresponding 
thermodynamic properties at the stagnation point. The stagnation -point heat -transfer 
rate is calculated by using the relation (ref. 12) 



where the constants rj . are given in the following table: 


(25) 


Species i 

^i 

co 2 

4.3102 X 10“ 4 

N 2 

3.6285 

02 

4.3102 

Ar 

5.4788 

He 

2.5038 

h 2 

1.2786 


The values of 77 * for C0 2 , N 2 , Ar, and H 2 are taken directly from table I of reference 12, 
whereas the value for He is obtained from table II of reference 13. Since a value of 
for 0 2 is not presented in reference 12, the r , ^ for 0 2 is set equal to that for C0 2 in the 
present study (ref. 13). 
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For case (3) (both free-stream and post-normal -shock flow conditions frozen), 
yg = yg = Hence, the perfect -gas normal-shock relations of reference 11 are used 

to determine the properties for static and for stagnation post -normal -shock conditions. 


Determining Quiescent -Acceleration-Gas Pressure 

An important parameter in the operation of an expansion tube is the quiescent pres- 
sure of the acceleration gas p^g. For a given condition in region (A) , it is p^g and 
acceleration -gas temperature and molecular weight which determine the extent of the 
unsteady flow expansion (that is, the velocity obtained in region (IT)); conversely, a par- 
ticular value of p^Q is required to obtain a given value of U5. In most studies in an 
expansion tube, U5 is one of the principal flow quantities. Therefore, to generate a 
desired U5, the required pjg must be known. 

The momentum equation for the normal shock moving into region (10) is equa- 
tion (2) where the subscript 10 replaces 1 and the subscript 20 replaces 2. At the 
interface of the acceleration gas and the test gas, it is required that P2Q equal pg and 
•U20 equal U5. Two cases providing a range of p^g are considered. The first is 
where the incident shock into the acceleration gas is traveling at the same velocity as the 
interface of the acceleration gas and the test gas (ref. 14). For this case, equation (2) in 
conjunction with the equation of state (where Zjq is unity, because of relatively low 
values of p^g, and T^q is equal to T]J becomes 


010 = 


— Ur 2 
T o 


The second case considered is where (ref. 5) 

y 10 + 1 

U S ,10 - -T— U 5 (27) 

Now, in the Langley 6 -inch expansion tube, helium is generally employed as the accelera- 
tion gas. For values of U s iq less than 10 km/sec or so, helium behaves as a perfect 
gas (ref. 15); hence, the perfect-gas relation of reference 11 gives, for y^g = 5/3, 

4 P 5 

P^O - 2 (28) 

AU Riv/r * _ 1 


where 


W 10 U s,10 

y m RT i 
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At U g jq greater than 10 km/sec, the helium acceleration gas departs from perfect- 

gas behavior because of ionization, and equation (28) is not valid. For example, at a 

M g 1Q of 20 , equation (28) yields a value of p 1Q that is approximately 1.2 times the 

value computed for real helium at a p^ of 1 N/m2. The effects of a real-helium 

acceleration gas are not included in the present program; hence, for U in greater 

s 3 lu 

than 10 km/sec, reference 15 should be consulted for more accurate values of p^ than 
those calculated by equation (28). 

DISCUSSION 

Flow quantities in region (a) are obtained by using basic measured inputs in the 
following combinations: 


Case ( 1 ) 

Pjj Tj, X i? and U g ^ 

Case ( 2 ) 

Pp Tj, X i? and P 2 

Case (3) 

Pl> Tjj P 4 , and X^ 


For each combination, the option of a standing shock at the secondary diaphragm may be 
exercised. The values in regions (T) and (Q) for case (1) where U g j is an input are 
compared, for several mixtures, to those calculated by using the program of reference 16. 
This comparison shows excellent agreement (all parameters within 0.5 percent) as 
expected, since both reference 16 and the present study employ references 6 and 7 as the 
source of thermodynamic properties. For case (2) where p 2 is an input, the values of 
P 2 calculated for case (1) where U Sj i is an input are in turn used as inputs; this cross- 
check shows excellent agreement between the results. For case (3) where p^ and T^ 
are inputs, the values in region ( 2 ) are observed to be in good agreement (U s j within 
0.2 percent) with those of reference 1 for helium driver gas when imperfect -gas effects 
are negligible. 

Combinations of measured input and options for obtaining stagnation -point conditions 
in region ( 5 ), when it is assumed that the flow conditions in region (A) are known, are 
summarized in the following table: 


— 

Case 

Measured 

input 

Unsteady 

expansion 

Postnormal 

shock 

(1) 

U 5 or p 5 

Equilibrium 

Equilibrium 

(2) 

U 5 or P 5 

Frozen 

Equilibrium 

(3) 

U 5 or P5 

Frozen 

Frozen 
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The flow quantities in region (?) for case (1) are compared with those obtained by 
using the real-air charts of reference 17. For this comparison, a 10-species (e*, Ar, 

N, N + , N2> O, 0 + , O2, NO, and NO + ) air model is employed in the present program, the 
air composition by volume being 78.08 percent N2, 20.95 percent O2, and 0.97 percent Ar. 
This comparison shows good agreement between the results obtained from the charts of 
reference 17 and from the present program. For the case where Pg is an input, the 
values of Pg calculated for the case where U5 is an input are in turn used as inputs. 
Results of this cross-check show good agreement. Quantities in regions (5^s) and (?^t) 
for thermochemical equilibrium are checked by comparing results with those obtained by 
using the program of reference 18. 

The large number of gas mixtures that may be used as test media coupled to the 
wide range of flow conditions that may be generated with these mixtures in the expansion 
tube results in an extremely large number of possible flow states. For this reason the 
present program, which requires relatively long computer times (to be discussed subse- 
quently), has not yet been exercised to generate a family of plots illustrating facility per- 
formance for several gases or gas mixtures. Instead, detailed sample calculations are 
presented. 

In performing a test in an expansion tube, the first consideration is to calculate the 
theoretical flow quantities for the chosen mode of operation. From the experimenter’s 
viewpoint, this is necessary in order to obtain the approximate magnitudes of velocity and 
pressure in the various regions. Knowledge of these magnitudes is required in the prep- 
aration of facility instrumentation for the test. Also, if several gas mixtures are to be 
used in a study, it may be desirable to determine conditions required to yield the same 
values of certain quantities for all gases (for example, velocity and/or pressure in the 
test section, or region (?))• 

As an example, the expansion tube is assumed to be operated in a low-to-moderate 
performance mode. The driver gas is helium at ambient temperature (T4 = 300 K) and a 
pressure of 34.474 MN/m2. The acceleration gas is also helium at ambient temperature. 
Three test gases are considered, air (N2~02-Ar mixture), pure CO2, and pure Ar. The 
Pj is varied from 861.9 N/m2 to 6.895 kN/m2 for all gases, and the Tj is ambient. 

The desired free-stream velocity U5 is 4.267, 4.877, 5.486, and 6.096 km/sec. For the 
experimenter, a knowledge is required of the approximate magnitudes of U s j and P2 
in region (?); Pg, Pg and q^ . in region (?); and the p^g to produce the desired 
U5. As shown in reference 18 a substantial reduction in computer time is realized by 
judicious selection of the number of species. Now, at such moderate performance, first 
and second ionization for air and for CO2 are negligible so that the species considered in 
the air model may be limited to Ar, N, N2, O, O2, and NO and the species considered for 
the CO2 are O, O2, C, CO, and CO2. For Ar, second ionization is assumed negligible and 
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the species considered are e", Ar, and Ar + . The inputs to the program and the procedure 
for varying the number of species are discussed in appendix C. The species inputs, in 
FORTRAN IV language, for the composition of these gases are 



Air 

O 

O 

CO 

Ar 

NUMSP 

6 

5 

3 

JINDX 

3 

2 

2 

IAR 

1 

0 

2 

IC 

0 

3 

0 

IC0 

0 

4 

0 

IC02 

0 

5 

0 

IN 

2 

0 

0 

m2 

3 

0 

0 

m0 

6 

0 

0 

10 

4 

1 

0 

102 

5 

2 

0 

PC02 

0.0 

1.0 

0.0 

PN2 

0.78 

0.0 

0.0 

P02 

0.21 

0.0 

0.0 

PAR 

0.01 

0.0 

1.0 

PMIX 

1.0 

1.0 

1.0 


The composition for the air model is Xj^ = 0.78, Xq^ = 0.21, and X^ r = 0.01. It 
should be noted that PMIX must satisfy the relation 

PMIX = PC02 + PN2 + P02 + PAR 

Symbols denoting position in species array (that is, IAR, IC, etc.) that are zero need not 
be included in the input. 

For a given p^, which in turn corresponds to a given set of flow quantities in 
region (A), four values of U5 are considered in this example; that is, four expansions 
from the given conditions in region (a) are to be performed. Hence, the inputs (appen- 
dix C) for this case are IREP = 1, U5I = 610, and NVEL = 4. For this example, the exis- 
tence of a standing shock at the secondary diaphragm is assumed (LD = 1) and only a 
thermochemical-equilibrium expansion from region (A) to region (IT) is considered 
(LG = 0). 

Pertinent quantities calculated within the present program for the three test gases 
considered in this example are given in tables I to III. 

The results for CO2 (table II) and for Ar (table III) were obtained on a Control Data 
6600 series computer, whereas the results for air (table I) were obtained on the slower 
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Control Data 6400 series computer. Three cases were also run for a 15-species CO2 
model and the same conditions as in table II for p^ = 6.895 kN/m 2 . These cases were 
run on a Control Data 6600 series computer and the computer times required per case 
are compared in the following table: 


Gas 

No. of species 

Computer time 
per case, 
sec 

Ar 

3 

24 

co 2 

5 

58 

C0 2 

15 

215 


At the moderate velocities in regions (2) and (T) considered, the computed flow quantities 
for the 5-species CO2 model and the 15-species CO2 model are the same. Therefore, 
species should be selected judiciously to minimize computer time without sacrificing 
accuracy in calculated flow quantities. 

Cases are also run, at an arbitrary p^ of 3.447 kN/m 2 , for gas mixtures repre- 
senting atmospheric models of Mars ^XQQg = 0.85 and X^ r = 0.15j, Venus = 0.95 

and X^2 = 0.05^, and Jupiter (X He = Xjjg = 0.5j. For the Mars model, the driver condi- 
tions are the same as in tables I to III and U5 is 5.486 km/sec. However, for the Venus 
model, maximum performance for resistance -heated helium driver gas in the Langley 
6 -inch expansion tube is utilized. This corresponds to p^ = 68.95 MN/m 2 and 
T4 = 600 K. The U5 for this Venus model is 10.97 km/sec. For the Jovian model, 
maximum performance with resistance -heated hydrogen driver gas (p^ = 68.95 MN/m 2 
and T4 = 600 is utilized and the U5 is 12.19 km/sec. Sample data printouts for 
these cases are presented in appendix C. 

For an actual test in an expansion tube, a comparison of the measured flow quanti- 
ties to the calculated flow quantities is desired. Such a comparison is illustrated in fig- 
ure 2 where preliminary experimental results for several tests in the Langley 6 -inch 
expansion tube are shown along with calculations from the present program. 

For these tests, the test gas is air and the driver gas is unheated helium at a P4 
of approximately 34.474 MN/m^. For a p^ of 3.447 kN/m2, the p^Q is 37 N/m 2 of 
helium and, for a p^ of 1.379 kN/m2, the p^Q is 3.3 N/m 2 of argon. In figure 2(a), 
the ratio is plotted against U g j and the solid line represents the calculated 

results shown in table I. From figure 2(a), it is observed that the U s \ inferred from 
the response of ionization gages located in the intermediate section is actually greater 
than the predicted U s j, being approximately 1.2 times the predicted values for both 
values of p^. This trend is consistently observed in the Langley 6 -inch expansion tube 
and is believed to result from a combination of several effects. Among these effects are 
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the following: (1) The cross-sectional area of the driver section is somewhat larger than 

the intermediate section. The data of figure 2(a) are obtained for a ratio of driver-section 
cross-sectional area to intermediate -section cross-sectional area equal to 1.173, whereas 
this ratio is assumed to be unity in the present program. As shown in reference 19, an 
increase in U s ^ is to be expected for an increase in this ratio. (2) The double - 
diaphragm mode of operation resembles that of a buffered -shock tube. Upon rupture of 
the upstream diaphragm, a shock wave propagates through the short chamber between the 
diaphragms. This shock wave heats the gas, contained within this chamber, just prior to 
rupture of the second diaphragm. Gains in U s ;[ obtained with buffered-shock tubes are 
illustrated in reference 19. (3) The relatively long opening times of the steel diaphragms 
may result in coalescing compressions and accelerating shock waves. It has been shown 
previously (ref. 20) that the finite opening time of the diaphragm produces faster shock 
waves for a given P4/P1 than simple shock-tube theory predicts. 

The Pg are plotted against U5 in figures 2(b) and 2(d) for Pj of 3.447 kN/m^ 
and 1.379 kN/m2, respectively; the Pg ^ are plotted against U5 in figures 2(c) and 2(e) 
for p^ of 3.447 kN/m^ and 1.379 kN/m2, respectively. The solid lines in figures 2(b) 
and 2(c) represent calculated results shown in table I. 

In figures 2(b) and 2(c), the experimental data show best agreement with the results 
calculated by the method which uses the measured U s j as an input, includes a standing 
shock at the secondary diaphragm, and assumes thermochemical equilibrium during the 
unsteady -expansion process. At the experimental values of p^ (fig. 2(b)) and p^ ^ 

(fig. 2(c)), the calculated velocity from table I is approximately 20 percent less than the 
measured velocity. This low prediction is expected, because of the difference in U s 1 
observed between theory and experiment in figure 2(a). The lower-density data of fig- 
ures 2(d) and 2(e) indicate a departure from an equilibrium expansion; for example, the 
calculated free-stream static pressure for an equilibrium expansion (fig. 2(d)) is approxi- 
mately 2.7 times the measured value, but this measured value is about 26 times as high 
as the calculated value for a frozen expansion (if it is assumed that a standing shock exists 
at the secondary diaphragm). 


CONCLUDING REMARKS 

A computer program written in FORTRAN IV language is presented which determines 
expansion -tube flow quantities for real test gases CO2, N2, O2, Ar, He, and H2, or mix- 
tures of these gases, in thermochemical equilibrium. This program permits, as input 
data, a number of possible combinations of flow quantities, which are generally measured 
during an expansion -tube test. The versatility of the program is enhanced by the inclusion 
of such effects as imperfect driver gas and shock reflection at the secondary diaphragm 
and by consideration of both a thermochemical -equilibrium flow expansion and a frozen 
flow expansion. 
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The usage of the program in preparing the expansion tube for testing with several 
gases is illustrated by sample calculations. Results from the program are also com- 
pared with preliminary data obtained experimentally in the Langley 6 -inch expansion tube. 
This comparison shows that the measured incident -shock velocity U Sj i and the mea- 
sured interface velocity U5 are approximately 1.2 times the calculated values, when 
driver-gas pressure and temperature are inputs. This discrepancy is attributed to 
several effects believed to be present in the Langley 6-inch expansion tube. Agreement 
between calculated and measured test -section flow quantities is significantly improved 
when measured incident -shock velocity is employed as an input. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 2, 1972. 
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APPENDIX A 


THERMODYNAMIC AND TRANSPORT PROPERTIES FOR 
ARBITRARY GAS MIXTURES 


Thermodynamic properties for real gases CO 2 , N 2 , O 2 , Ar, He, and H 2 , or mix- 
tures of these gases, are obtained by the procedure of references 6 and 7 and include 
dissociation and first and second ionization. Basic assumptions used in obtaining these 
properties are 

(1) The mixture is composed of ideal gases 

(2) For diatomic species the rigid-rotor harmonic-oscillator model is used with 

vibrational -rotational corrections for each electronic configuration 

(3) Only electronic levels with principal quantum number less than or equal to 5 

are included 

The procedure of references 6 and 7 is based upon the free -energy -minimization 
method of reference 21. For a given pressure and temperature, the free energies for 
individual species are computed from the partition function of statistical mechanics. The 
equilibrium composition is then obtained by minimization of the free energy. When the 
composition for a given pressure and temperature is determined, the corresponding ther- 
modynamic properties are computed directly. 

The 5 components and 26 species considered in references 6 and 7 are listed in 
appendix C. These 26 species permit modeling of Earth, Mars, and Venus atmospheres 
but not of proposed Jovian atmospheres which are believed to consist primarily of helium 
and hydrogen. Therefore, the components H and He and the species H, H 2 , H + , He, He + , 
and He ++ are also included in the present study. However, no provision is included for 
combinations of hydrogen with any of the other species. The energy -level constants for 
the 26 species considered in references 6 and 7 are tabulated in reference 6, whereas 
the energy-level constants for the 6 additional species are obtained from reference 22. 


The procedure of references 6 and 7 is included in the present study as a subroutine 
and is referred to herein as ROGO. The basic inputs to ROGO, other than the energy - 
level constants, are pressure, temperature, and initial estimates of the species concen- 
trations. (In the free -energy-minimization method for determining the composition of a 
gas mixture, the initial estimates of Yi are somewhat arbitrary.) In the present study, 
the partial pressures of the quiescent test gas in the intermediate chamber are assumed 
known. The corresponding mole fractions for these gases are obtained from Dalton’s law 


Xi = 


£i 

Pmix 


(Al) 
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APPENDIX A - Continued 


and the molecular weight is given by 


w o = y XiWi 
i 


(A2) 


The initial estimates for the species concentrations for CO 2 , N 2 , O 2 , Ar, He, and H 2 are 
obtained from 


Yi = 


1_ W 


(A3) 


For the remaining 26 species the initial estimates for the concentrations are set equal 
to 10 - 20. (if for O 2 is not known, the initial estimate is set equal to 10"5.) Ther- 
modynamic properties obtained from ROGO, for a given p and T, are a, h, sW 0/ /R, 

x i> z * > rE’ and p ‘ 


Combinations of input thermodynamic properties other than p and T are required 
in the present study. An iterative -interpolation scheme is derived in reference 18 so as 
to enable determination of thermodynamic properties from combinations of h, p, sW 0/ /R, 
and p. This scheme is referred to herein as FIND (I), and the inputs for a given I are 
as follows: 


I 

Inputs to FIND (I) 

1 

P, P 

2 

p, h 

3 

p, sW 0 /R 

4 

h, p 


A detailed discussion of FIND (I) is presented in reference 18. 

The viscosity p is calculated by using the expressions presented in reference 18. 
Since He and H 2 are also included in the present study, expressions for the viscosity of 
these species are derived. The viscosity for He is obtained by applying a curve fit to the 
results of references 23 to 25 for temperatures to 5000 K. The resulting expression is 

p = 3.92 x 10 _7 T 0 * 7 (A4) 

The expression for Ar is (ref. 18) 

p = 3.33 X io- 7 T 0 * 739 (A 5) 

Simple expressions for the viscosity, in the form 

p = (b 0 + bjT + b 2 T 2 ) X 10" 7 (A6) 

are obtained for the diatomic species N 2 and O 2 and for the triatomic species CO 2 in ref- 
erence 18. For convenience, the constants of equation (A6) for these species along with 
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those for H2 obtained from the results of references 25 to 27 are presented in the fol- 
lowing table: 


Gas 

T, K 

b 0 

b r 

b 2 

co 2 

150 to 1000 
1000 to 5000 

-9.7303 x 10' 1 
1.5029 X 10 2 

5.5222 x 10 _1 
2.4890 X 10" 1 

-1.6139 x 10" 4 
-6.1747 X 10' 6 

N 2 

150 to 1000 
1000 to 5000 

2.2992 X 10 1 
1.6388 X 10 2 

5.5586 x 10 _1 
2.4601 x 10" 1 

-1.8436 X 10 -4 
-7.7388 X 10' 6 

°2 

150 to 1000 
1000 to 5000 

1.9939 x 10 1 
2.3030 x 10 2 

6.6236 x 10 _1 
2.5597 X 10 -1 

-2.0410 X 10‘ 4 
-2.2643 X 10" 6 

H 2 

150 to 1000 
1000 to 5000 

2.3877 x 10 1 
6.3197 X10 1 

2.1918 X 10 _1 
1.4093 x 10" 1 

-4.4304 X 10~ 5 
-8.6387 X IQ' 6 


These expressions should provide reasonably accurate (within 10 percent or so) values 
of ju for T less than 5000 K, with the accuracy decreasing with further increase in T. 
For T greater than 8000 K, \± is not calculated in the present program. 
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APPENDIX B 


THERMODYNAMIC PROPERTIES FOR IMPERFECT 
HELIUM AND HYDROGEN 

The driver gas for the Langley 6 -inch expansion tube is usually helium or hydrogen. 
For unheated and resistance -heated driver gas, pressures up to 70 MN/m 2 and tempera- 
tures from 300 K to 600 K are obtained. For such conditions, the assumption of perfect 
helium or hydrogen is inadequate because of intermolecular -force effects (for example, 
see refs. 28 and 29). The relations used to calculate imperfect-gas effects at these 
driver -gas conditions are taken from reference 28 and are based on the virial form of 
the equation of state which, where interactions involving four or more molecules are 
neglected, can be written as 

p = p^-TZ = p J^-T[l + pB(T) + p 2 C(T)] (Bl) 

For helium, the virial coefficients are given by 

B(T) = 3.3565 x 10' 6 (15.8922 - In T) 3 

- 2.0085 x 10' 3 exp[-(3.7156 X 10 _3 )t| (B2) 

and 

C(T) = 5.6330 x 10" 12 (15.8922 - In T) 6 (B3) 

These expressions for B(T) and C(T) are employed in the imperfect-helium study of 
reference 28 for a temperature range of 200 K to 15 000 K. 

For hydrogen a curve fit is applied to the second-virial -coefficient results of ref- 
erence 30. These results include a temperature range of 273 K to 473 K and the curve - 
fit expression is 

B(T) = 1.6994 x 10' 3 T 1//4 (B4) 

Although equation (B4) is obtained for temperatures to only 473 K, this expression is 
assumed to be valid for temperatures to 600 K. (This temperature represents the maxi- 
mum value of T4 considered herein.) To obtain an expression for C(T) for hydrogen, 
the results of reference 29 are used. In reference 29 compressibility factors are calcu- 
lated by using the Beattie -Bridgeman equation of state for temperatures from 30 K to 
1000 K and for pressures from 10 kN/m 2 to 10 2 MN/m 2 . Now, for pressures to 
10 MN/m 2 , the compressibility factors obtained from equations (Bl) and (B4), where 
C(T) = 0, are within 1 percent of the values presented in reference 29 for temperatures 
from 300 K to 1000 K. However, at higher pressures (lO MN/m 2 < p < 10 2 MN/m 2 ), these 
compressibility factors are as much as 20 percent less than those of reference 29, since 
interactions involving three molecules are ignored (that is, C(T) = 0), Hence, values of 
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C(T) are calculated which bring the compressibility factors from equation (Bl) into agree- 
ment with those of reference 29 for a pressure of 10 2 MN/m 2 . For a temperature range 
of 300 K to 1000 K, these values of C(T) vary from 2.002 x 10“ 4 to 2.195 x 10 -4 . With no 
appreciable loss of accuracy, C(T) is set equal to the constant 2.1 x 10 -4 for the tem- 
perature range of 300 K to 600 K. The B(T) of equation (B4) and this constant value of 
C(T) result in compressibility factors within approximately 1 percent of those of refer- 
ence 29 for the temperature range of 300 K to 1000 K and the pressure range of 10 2 kN/m 2 
to 10 2 MN/m 2 . 

The thermodynamic relations, in terms of the virial coefficients, for hW 0 /RT, 
sW 0 /R, c v W 0 /R, and a 2 are (ref. 28) 


^ JhW\ 


RT 


RT / 


perf 


+ P 


B(T) - T 

dT 


2C(T) - T 


dC(T) 


SW, 


o _ ( c v Wq 
R 


perf 


In T - In p - p 


B(T) + T 


dT 


fl!| 

2 


dT 


C(T) + T ^2111 
.dT 


sW ( 

B /ref 


(B5) 


(B6) 


c v W 0 


c v W 0 


perf 


- T<p 


2 dB(Tl , T d^B(T) 

_ dT dT 2 J 




2 d.C(T) , T d z C(Tj 
dT dT 2 


(B7) 


,2 = 


8 p/t c v p 2 Wp 


where, from differentiation of equation (Bl), 



B(T) + T 


dB(T) 

dT 


+ P‘ 


C(T) + T 

dT J / 


(|)t = |r[l + 2 PB(T) + 3p 2 C(T)] 


(B8) 


(B9) 


(B10) 


The ( sW o/ R )ref e ^ uation (B6) I s obtained for p = 1.01325 x 10 2 kN/m 2 and 
T = 300 K by using reference 10 for both helium and hydrogen. For helium, a value of 
( sW o/ R ) re f = 4 - 8024 is calculated and, for hydrogen, ( sW o/ R ) re f = -1«0363. (it should 
be noted that ( s W 0 /R) re j is omitted in reference 28; hence, the values of sW 0 /R 
plotted for helium in reference 28 differ from equation (B6) by a constant ( sW o/ R ) re f^ 
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In the thermodynamic relations shown in equations (B5) to (BIO), the independent 
variables are T and p. Hence, in generating an array of thermodynamic properties in 
region (3) for a range of T3 and a known S3W3/R (see section entitled "Determining 
Flow Quantities in Region the corresponding Pg must be known. These Pg are 

determined iteratively from the implicit relation given by equation (B6). With Tg known 
and the corresponding Pg determined, corresponding values of Z and p, h, and a 
may be calculated from equations (Bl), (B5), and (B8), respectively. 
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APPENDIX C 


COMPUTER -PROGRAM INPUTS, LISTING, AND FLOW CHART 
WITH SAMPLE DATA PRINTOUTS 


The present program is written in FORTRAN IV language for Control Data 6000 
series computers. Minimum machine requirements are 70 000 octal locations of core 
storage. The FORTRAN NAMELIST capability is used for data input with INP as the 
NAMELIST name. The units for the inputs which are physical quantities are given in the 
section entitled "Symbols." The program symbols and a brief description of the inputs 
necessary to utilize the computer program are listed as follows: 

Program symbol 

Description 

PI 

Pressure of quiescent test gas in region (T) 

T1 

Temperature of quiescent test gas in region (T) 

US1 

Incident -shock velocity into region (T) 

P2 

Static pressure in region (IT) 

P4 

Driver-gas pressure in region (T) 

T4 

Driver-gas temperature in region (T) 

U5 

Velocity in region (IT) 

P5 

Static pressure in region (IT) 

PC02 • 

Partial pressure of CC >2 

PN2 

Partial pressure of N 2 ' 

P02 

Partial pressure of O 2 

PAR 

Partial pressure of Ar 

PHE 

Partial pressure of He 

PH2 

Partial pressure of H 2 

PMIX 

Mixture pressure 

TW 

Model surface temperature 
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Program symbol 
BNR 
RUN 
NUMSP 
JINDX 
IAR 
IC 
IC0 
IC02 
IN 
IN 2 
IN0 
10 
102 
IHE 
IH 
IH2 

NDRIV 

LB 


Description 

Model nose radius 
Facility test number 

Number of species considered (30 maximum) 

Number of components considered (8 maximum) 

Position in species array of Ar 
Position in species array of C 
Position in species array of CO 
Position in species array of CO2 
Position in species array of N 
Position in species array of N2 
Position in species array of NO 
Position in species array of O 
Position in species array of O2 
Position in species array of He 
Position in species array of H 
Position in species array of H2 
NDRIV = 0 denotes helium driver gas 
NDRIV = 1 denotes hydrogen driver gas 

LB = 0 denotes inputs p^, Tj_, and U s j used to find region (2) 
quantities 

LB = 1 denotes inputs p^, Tj, and Pg used to find region (IT) 
quantities 
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LB = 2 denotes inputs p^, Tj, p^, and T4 used to find 
region (T) quantities 
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Program symbol 
ISTET 


LF 


LD 


LG 


IREP 


U5I 

NVEL 

NUMUS 


Description 

ISTET = 0 denotes only quantities in regions ( 2 ) and (2^s) 
determined 

ISTET = 1 denotes all expansion -tube flow quantities determined 

LF = 1 denotes U5 is basic input in region (T) 

LF = 2 denotes Pg is basic input in region (5) 

LD = 0 denotes no standing shock at secondary diaphragm 

LD = 1 denotes existence of standing shock at secondary 
diaphragm 

LD = 2 denotes both cases (LD = 0 and LD = 1) are performed 

LG = 0 denotes frozen expansion is not performed 

LG = 1 denotes frozen expansion is performed 

IREP = 0 denotes only a single value of U5 is of interest for 
given region (A) quantities 

IREP = 1 denotes several U5 of interest for given region (A) 
quantities 

Velocity increment for IREP = 1 

Total number of U5 of interest for IREP = 1 (10 maximum) 

Number of U s j in array for LB = 2 


To minimize the number of inputs required for running cases on the computer, 
inputs are assigned values within the program as follows: 


Program symbol 


Assigned value 


TW 

300 

BNR 

0.0127 

RUN 

1.0 

NDRIV 

0 
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Program symbol 
ISTET 
LD 
LG 
IREP 
U5I 
NVEL 
NUMUS 


Assigned value 
1 
2 
1 
0 
0 
0 

16 


These values may be changed from their assigned values by inclusion in the NAMELIST 
INP. For a given LB, only the basic parameters Pp Tp and U s ^ (LB = 0), p 2 
(LB=l),or p^ and T 4 (LB = 2) need be included in INP. Similarly, for a given LF, 
only U 5 (LF = 1) or p^ (LF = 2) need be included in INP. Only the partial pressures 
of the quiescent test gas in region (T) need be included in conjunction with PMIX. It 
should be noted that the condition 

PMIX = PC02 + PN2 + P02 + PAR + PHE + PH2 

must be satisfied. Along with LB and LF, NUMSP and JINDX must also be specified for 
each case. However, only the positions in the species array (for example, IAR, IC, IC0, 
and so forth) of the species being considered need be included. 

The 26 species and 5 components considered in references 6 and 7 and employed 


in the present program are 

as follows 




Species 


Components 

e" 

O 

c 

e" 

Ar 

o + 

,c + 

Ar 

Ar + 

0 ++ 

c ++ 

N 

Ar ++ 

o' 

C" 

O 

N 

°2 

CO 

C 

N + 

° 2 + 

co + 


N ++ 

V 

CN 


n 2 

NO 

co 2 


n 2 + 

Z 

O 

+ 
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APPENDIX C - Continued 


The thermodynamic data for these species, which are tabulated in reference 6, are read 
into the computer program from cards. A listing of these cards is presented in refer- 
ence 18. An additional 6 species and 2 components considered in the present study are 
as follows: 



Species 

Components 

H 

He 

He 

H + 

He + 

H 

h 2 

He ++ 



The thermodynamic data for these species are obtained from reference 22 and are also 
read into the computer program from cards. 

A listing of this program, including subroutines, comments, sample input, and 
thermodynamic data for a mixture of He and H 2 , is reproduced in the following pages. 
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print 75. W0.XC02.XN2.X02.XAR.XHE.XH2 

2MIX=1 «E+4*(7.2906*XAR+10.21 09*XC02+7 .72 1 1 *XN2+7 .4243*X02+ 1 .5988*X 
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IF (P2( 1 >«GT.P2C2> ) 
IF ( P3 ( 1 >.GT*P3<2> ) 
P31 *P3 ( 1 ) 

P32=P3(M) 

P21 =P2 ( 1 ) 
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APPENDIX C — Continued 


The basic subroutines of this program are as follows: 

(1) ROGO - computes thermodynamic quantities a, h, sW 0 /R, Xi, Z*, yg, 

and p for given p, T, and W 0 

(2) FIND - iterative -interpolation procedure for obtaining thermodynamic quanti- 

ties from ROGO with combinations: 

1 . p and p 

2. p and h 

3. p and sW 0 /R 

4. h and p 

(3) VISC - computes p for given T and X^ 

(4) BDT - computes B and C for helium or hydrogen for given T 

(5) SOLUT - given (p2>U2) and (P3JU3), finds solution to curves 

(6) SC - iterative procedure for solving conservation relations for a moving 

normal shock 

(7) SNS - iterative procedure for solving conservation relations for a standing 

normal shock 


(8) SIMR - computes 



by Simpson's rule 


A flow chart of this program is given in the following pages. 
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APPENDIX C - Continued 



HELIUM 


HYDROGEN 


CONSTANTS CONSTANTS 


CALL 


PRINT 

HEADINGS 


CALL 
SI MR 


LB = 2? 


CALL 


CALL SC 


CALL 

SOLUT 


CALL 


CALL 

ROGO 


T 2 > 8000? 


NDRIV 
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APPENDIX C - Continued 


The procedure for running a program with all 26 species for CO 2 , N 2 , 0 2 , and Ar 
is discussed in detail in reference 18. Also, modifications required when the number of 
species is varied are illustrated by example in reference 18. For convenience, an addi- 
tional example is now presented. A mixture of CO 2 and Ar, which may represent the 
Mars atmospheric model, is considered with second ionization neglected. Now, the num- 
ber of atoms of a particular component per particle of a species A(I,J) must be desig- 
nated. For electrons, a charge constraint is used instead of a mass constraint. Con- 
structing a table of A(I,J) for this C02~Ar mixture gives 


12 3 4 



^^^^Component (J) 
Species (I) 

e " 

Ar 

0 

C 

1 

e" 

1 

0 

0 

0 

2 

Ar 

0 

1 

0 

0 

3 

Ar" 1 " 

-1 

1 

0 

0 

4 

O 

0 

0 

1 

0 

5 

o + 

-1 

0 

1 * 

0 

6 

0 " 

1 

0 

1 

0 

7 

°2 

0 

0 

2 

0 

8 

<V 

-1 

0 

2 

0 

9 

°2~ 

1 

0 

2 

0 

10 

C 

0 

0 

0 

1 

11 

C + 

-1 

0 

0 

1 

12 

C" 

1 

0 

0 

1 

13 

CO 

0 

0 

1 

1 

14 

CO+ 

-1 

0 

1 

1 

15 

co 2 

0 

0 

2 

1 


The sequence 0, 1, 0, 0 for species Ar denotes that the only component contained in this 
species is Ar. (A zero denotes that a particular component does not appear in that 
species.) In the sequence -1, 1, 0, 0 for Ar + , the -1 denotes* that Ar + is a positive ion, 
and so forth. For this mixture, NUMSP = 15, JINDX = 4, IAR = 2, IC = 10, IC0 = 13, 
IC02 = 15, 10 = 4, and 102 = 7. 

Only the species decks for the 15 species considered are included, and these are 


placed in the order that they appear in the 
illustrated by the following listing: 

A(I,J) 

table. 

The 

A(I,J) 

are read in as 

1 

0 

0 

0 

0 

1 

• 0 

0 

-1 

1 

0 

0 0 

0 

1 

0 

-1 

0 

1 

0 

1 

0 

1 

0 

0 

0 

2 0 

-1 

0 

2 

0 

1 

0 

2 

0 

0 

0 

0 

1 

-1 

0 0 

1 

1 

0 

0 

1 

0 

0 

1 

1 

-1 

0 

1 

1 

0 0 

2 

1 
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APPENDIX C - Continued 


The first four numbers represent the A(I,J) sequence, from the A(I,J) table, for e”; 
the next four numbers represent the A(I,J) sequence for Ar; and so forth. The sequence 
of cards following the main program and subroutine is as follows: 



Sample data printouts are included that illustrate program usage for hypothetical 
tests in an expansion tube with a Mars model atmosphere at moderate velocity and with 
Venus and Jovian model atmospheres at high velocity. 

The headings for the various flow regions correspond to those in the section 
entitled "SYMBOLS.” The units of these flow quantities are as given in the definitions 
of the symbols, except where explicitly specified in the printout. The sample printouts 
are as follows: 
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P PHO T H PS1 US1 

3Q3E+01 9.847E+0C 7.316E+00 6.942E-01 7.815E+00 2.146E»03 
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TABLE I.- CALCULATED EXPANSION -TUBE FLOW QUANTITIES FOR SIX-SPECIES 
AIR MODEL IN THERMOCHEMICAL EQUILIBRIUM 
[j8% N 2 , 21% 0 2 , 1% Ar, by volume; P4 = 34.474 MN/m 2 ; 

T 4 = T 4 = T 10 = 300 K; W 4 = W 10 = 4.003 kg/kmol; 

W 4 = 28.97 kg/kmoTJ 


U 5’ 

km/sec 

Ps ’ 0 
kN/m 2 

t 5 ,k 

Z*5 

V E,5 

m 5 

N Re,5> 

m _ l 

P5 t> 
kN/m 2 

^5,t> 

MW/m 2 

PlO’ 

N/m 2 

P x = 861.9 N/m 

2 ; p 2 = 

58.25 kN/m2 ; 

T 2 = 2766 K; Z* 2 = 

1.013; 

U S) 1 = 2.579 km/sec; 


u 2 = 2.223 km/sec; P 2 s = 326.2 kN/m2 ; T 2 g = 3680 K; 

Z* 2 ,s = 1-OT1 

4.267 

3.727 

1924 

1.000 

1.271 

5.094 

4.547 X10 5 

119.2 

12.74 

120.2 to 123.3 

4.877 

1.396 

1545 

1.000 

1.301 

6.420 

2.803 

72.62 

12.03 

31.20 to 35.65 

5.486 

.4560 

1184 

1.000 

1.323 

8.183 

1.587 

39.18 

10.64 

7.476 to 9.248 

6.096 , 

.1229 

852.1 

1.000 

1.349 

10.61 

.8080 

18.12 

8.663 

2.016 to 2.028 

p x = 3.447 kN/m 2 ; p 2 = 

190.4 kN/m 2 ; 

T 2 = 2495 K; Z* 2 = 

= 1.002; 

U g ^ = 2.345 km/sec; 


U 2 = 1.991 km/sec; P 2 s = 940.4 kN/m 2 ; T 2 g = 3450 K; 

z *2,s * i- 028 

4.267 

5.871 

1305 

1.000 

1.315 

6.079 

13.59 X10 5 

276.8 

18.03 

189,4 to 194.2 

4.877 

1.706 

962.5 

1.000 

1.339 

8.019 

7.316 

142.5 

15.97 

38.11 to 43.55 

5.486 

.3869 

652.9 

1.000 

1.370 

10.83 

3.536 

60.29 

12.71 

6.343 to 7.847 

6.096 

.05869 

386.8 

1.000 

1.396 

15.48 

1.475 

19.06 

8.656 

.9621 to .9679 

Pj = 6.895 kN/m 2 ; p 2 

335.9 kN/m 2 ; T 2 = 

2291 K; Z* 2 

= 1.001; 

U c - = 2.207 km/sec; 

S,1 


U 2 = 1.861 km/sec; P 2 s = 1*559 MN/m 2 ; T 2 g = 

3243 K; 

Z*2,s = 

1.013 

4.267 

6.016 

1002 

1.000 

1.336 

6.885 

21.20 X10 5 

369.4 

20.11 

194.1 to 199.0 

4.877 

1.415 

688.3 

1.000 

1.366 

9.389 

10.52 

165.2 

16.79 

31.62 to 36.13 

5.486 

.2283 

416.1 

1.000 

1.394 

13.45 

4.541 

55.82 

12.04 

3.742 to 4.630 

6.096 

.01908 

205.1 

1.000 

1.401 

21.23 

1.481 

11.69 

6.712 

.3129 to .3148 
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TABLE H.- CALCULATED EXPANSION -TUBE FLOW QUANTITIES FOR FIVE -SPECIES 
PURE C0 2 IN THERMOCHEMICAL EQUILIBRIUM 
jp 4 = 34.474 MN/m2; T 4 = T x = T 10 = 300 K; W 4 = W 10 = 4.003 kg/kmol; 

W 4 = 44.011 kg/kmolj 


U 5 . 

km/sec 

Ps> 

kN/m 2 

t 5 ,k 

— 

r7* 

Z 5 

y E,5 

M 5 

N Re,5> 

m-1 

P5,t> 

kN/m 2 

45,t> 

MW/m 2 

P 10’ 

N/m 2 

Pi - 

861.9 N/m2 ; p 2 = 

= 79.10 kN/m 2 

T 2 = 2275 K; Z* 2 = 

1.032; U s i as 2.371 km/sec; 


U 2 = 2.169 km/sec; P 2 s = 797.9 kN/m 2 ; T 2 g = 2990 K; Z 

*2,s = 1 - 138 

4.267 

2.131 

1933 

1.019 

1.107 

6.649 

4.016 xlO 5 

101.2 

13.71 

68.75 to 70.51 

4.877 

.7220 

1749 

1.008 

1.122 

7.977 

1.864 

49.99 

11.70 

16.13 to 18.44 

5.486 

.2251 

1534 

1.002 

1.147 

9.506 

.8214 

22.62 

9.478 

3.689 to 4.565 

6.096 

.06330 

1291 

1.001 

1.169 

11.41 



.3428 

9.352 

7.272 

1.038 to 1.044 

Pi 3 

3.447 kN/m 2 ; p 2 = 253.2 kN/m 2 ; T 2 = 

2084 K; Z* 2 

= 1.009; 

U c i = 2.130 km/sec; 


u 2 = 

1.927 km/sec; 

P 2,s = 

2.230 MN/m 2 ; T 2 s = 

2844 K; 

Z* 2)S = 1-070 

4.267 

3.336 

1479 

1.001 

1.162 

7.485 

10.08 XlO 5 

210.7 

18.33 

107.6 to 110.4 

4.877 

.9050 

1226 

1.000 

1.173 

9.352 

4.267 

90.07 

14.80 

20.22 to 23.10 

5.486 

.2127 

985.6 

1.000 

1.182 

11.69 

1.608 

33.33 

11.00 

3.486 to 4.313 

6.096 

.04114 

759.3 

1.000 

1.197 

14.71 

.5373 

10.33 

7.390 

.6744 to .6784 

pi 31 

6.895 kN/m 2 ; p 2 

= 441.6 kN/m 2 ; T 2 = 

1929 K; Z* 2 ; 

= 1.003; 

U s j = 1.993 km/sec; 


U 2 = 1.793 km/sec; 

P 2 s = 3.634 MN/m 2 ; T 2 s = 

2713 K; ! 

z *2,s= I- 042 

4.267 

3.300 

1185 

1.000 

1.175 

8.321 

14.41 X10 5 

260.3 

19.61 

106.5 to 109.2 

4.877 

.7534 

946.4 

1.000 

1.184 

10.60 

5.447 

97.17 

14.96 

16.84 to 19.24 

5.486 

.1405 

722.3 

1.000 

1.200 

13.56 

1.800 

30.05 

10.23 

2.304 to 2.850 

6.096 

.01944 

511.0 

1.000 

1.227 

17.71 

.5133 

7.255 

6.105 

.3188 to .3206 
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TABLE m.- CALCULATED EXPANSION -TUBE FLOW QUANTITIES FOR 
THREE -SPECIES PURE Ar IN THERMOCHEMICAL EQUILIBRIUM 
^> 4 = 34.474 MN/m2 ; T 4 = Tj = T 10 = 300 K; W 4 = Wi 0 = 4.003 kg/kmol; 

Wj = 39.944 kg/kmolj 


U 5> p 5> 

km/sec N/m 2 


T 5 - k Z*5 



„ 10 ’ 
/m 2 MW/m 2 N/m 2 


Pj * 861.9 N/m 2 ; p 2 = 86.82 kN/m 2 ; T 2 = 7723 K; Z* 2 = 1.001; U s l = 2.893 km/sec; 
U 2 = 2.153 km/sec; P 2 s = 192-9 kN/m 2 ; T 2 s = 9779 K; Z* 2 s = 1.012 



2659 

1.000 

1.667 

4.443 

3.820 X10 5 

178.8 

1653 

1.000 

1.667 

6.440 

3.041 

114.4 

884.7 

1.000 

1.667 

9.903 

2.125 

56.70 

354.6 

1.000 

1.667 

17.38 

1.177 

17.76 



180.6 to 185.3 
38.13 to 43.57 
5.863 to 7.254 
.5961 to .5997 


Pl = 3.447 kN/m 2 ; p g * 275.4 kN/m 2 ; T 2 = 6267 K; Z* 2 = 1.000; U s l = 2.583 km/sec; 
U 2 = 1.907 km/sec; p 2 g = 518.2 kN/m 2 ; T 2s = 8065 K; Z* 2)S = 1.001 


4.267. 6245 1389 1.000 1.667 6.146 

4.877 1106 695.1 1.000 1.667 9.931 

5.486 76.52 238.8 1.000 1.667 19.06 

6.096 .1675 20.61 1.000 1.667 72.08 



381.4 

29.04 

176.3 

24.70 

44.94 

15.48 

1.407 

3.362 


201.4 to 206.6 
24.70 to 28.22 
1.255 to 1.552 


P x =• 6.895 kN/m 2 ; p 2 = 477.2 kN/m 2 ; T 2 = 5469 K; Z* 2 » 1.000; U s l = 2.405 km/sec; 
U 2 * 1.771 km/sec; p 2 g = 871.6 kN/m 2 ; T 2 g = 6999 K; Z* 2 s = 1.000 


4.267 4775 

4.877 476.8 


5.837 

.07358 


873.0 1.000 

347.3 1.000 

59.69 1.000 1.667 
10.37 1.000 1.667 



22.59 X10 5 
12.80 
3.770 
1.108 



31.14 154.0 to 158.0 

22.60 10.65 to 12.17 

8.498 .09569 to .1184 

3.141 .001206 to .001213 






































I nterface 



Distance ► 

(b) Reflected shock at secondary diaphragm. 
Figure 1.- Concluded. 
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(a) Incident-shock velocity as function of pressure ratio across primary diaphragm. 

Figure 2.- Comparison of calculated and measured expansion -tube flow quantities 
for unheated helium driver gas and air test gas. 



m 


Symbols denote measured quantities 


Calculated 

Standing shock; equilibrium expansion 

Equilibrium expansion 

Standing shock; equilibrium expansion 

Frozen expansion 

Standing shock; frozen expansion 


V T 


(b) Free -stream static pressure as function of free -stream velocity; 

Pj = 3.447 kN/m2. 

Figure 2.- Continued. 




(c) Pitot pressure as function of free-stream velocity; = 3.447 kN/m 

Figure 2.- Continued. 







(e) Pitot pressure as function of free-stream velocity; = 1.379 kN/m^, 

Figure 2.- Concluded. 
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